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Рассмотрена осесимметричная задача теории упругости о кручении недеформируемым круглым штампом 
упругого полупространства с неоднородным по глубине покрытием периодической структуры. Покрытие 
представляет собой многослойный пакет с чередующимися мягкими и жёсткими слоями, при этом границы 
слоёв могут быть как чёткими (кусочно-постоянное изменение модуля сдвига), так и сглаженными (непре- 
рывно-неоднородные покрытия). Построено приближённое аналитическое решение задачи высокой точно- 
сти, эффективное для любых толщин покрытий. Изучено влияние числа слоёв и характера неоднородности 
покрытия на трансформанту ядра интегрального уравнения и на распределение контактных напряжений под 
штампом. Показано, что при определённых параметрах задачи достигается существенное различие между 
результатами для непрерывного и скачкообразного характера изменения модуля сдвига по глубине. 
Ключевые слова: кручение, многослойные покрытия, слоистые композиты. 


Введение. Многослойные покрытия активно используются для создания современных жаропроч- 
ных, абразиво- и эрозионностойких покрытий на поверхности элементов газовых и паровых тур- 
бин, деталей машин. Разработан ряд технологий, позволяющих создавать слоистые композиции 
из чередующихся слоёв различных материалов, в которых толщина слоя может составлять менее 
100 нм при общем их количестве до нескольких тысяч [1]. Это газотермическое напыление (в 
частности, вакуумно-плазменное) и различные способы эпитаксиального нанесения (осаждения 
из газовой или водной среды). В качестве компонентов покрытия могут быть использованы раз- 
личные металлы и сплавы, а в случае эпитаксии — полимерные материалы. 

В работе [2] рассмотрены многослойные покрытия периодической структуры (чередующи- 
еся слои из алюминия и палладия) со сглаженными границами между слоями. На основании экс- 
периментальных результатов показана перспективность использования подобных покрытий для 
защиты от разрушения при механическом воздействии. 

В настоящей работе анализируются механические характеристики контактного взаимодей- 
ствия непрерывно-неоднородных и многослойных покрытий периодической структуры, возника- 
ющие при кручении упругого однородного полупространства с покрытием недеформируемым 
круглым штампом. 

Задача о кручении однородного упругого полупространства круглым штампом впервые 
была сформулирована и решена в динамической постановке Рейснером и Сагочи [3]. Снеддон [4], 
с использованием техники интегральных преобразований, свёл данную задачу к решению инте- 
грального уравнения. 


“Работа выполнена при поддержке Министерства образования и науки РФ (ГК № 11.519.11.3015, соглашение 
№ 14.В37.21.1131) и РФФИ (гранты № 12-07-00639-а, 13-07-00954-а). 
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В. Д. Грилицкий [5] построил решение задачи о кручении круглым штампом изотропной 
двухслойной среды и ортотропного упругого слоя в виде степенного ряда от отношения толщины 
первого слоя к радиусу штампа. 

В работе [6] построено решение задачи о кручении для упругого трансверсально- 
изотропного полупространства с неоднородным по глубине трансверсально-изотропным покрыти- 
ем. Решение эффективно для всего набора физических и геометрических параметров задачи. 
Постановка задачи и построение решения. Недеформируемый круглый штамп с плоским ос- 
нованием жёстко сцеплен с верхней гранью Г упругого неоднородного полупространства “2, со- 
стоящего из неоднородного по глубине покрытия толщины НЯ и однородного полупространства 
(подложки). С полупространством связана цилиндрическая система координат г, ф, 2 о6ь 2 
нормальна плоскости Г и проходит через центр штампа, координата г отнесена к радиусу штампа 
а. Штамп контактирует с полупространством по поверхности 2 =0, г<1. К штампу приложен 


крутящий момент //, ось которого совпадает с осью 2. Под действием этого момента штамп по- 
вернётся относительно оси 2 на угол =, вызвав деформацию кручения полупространства &. 
Модуль сдвига полупространства изменяется с глубиной по закону 


6°Р (2), -Н<2<0 
С°, —©<2<-Н 
где /(2) — непрерывно-дифференцируемая или кусочно-постоянная функция. 


Учитывая, что слой подвержен лишь деформации скручивания, ненулевым является лишь 
смещение вдоль оси ф: 


С(2)= ‚ ($ = СОП5Е, 


ПО: ИО Рай 


Вне штампа плоскость Г не нагружена. При сделанных предположениях граничные условия 
имеют вид: 
=0, />1, 


Е. ГЫ, 


т 
2=0: в =Т. =0, 3 ^ 
Ф 

При г > © И 2 -> —« напряжения исчезают. 
Считаем, что перемещения и напряжения сопрягаются на границе между покрытием 


(верхний индекс (1)) и подложкой (верхний индекс (2)): 


2=-Н: т®=то, и =и0. 


Требуется определить закон распределения контактных касательных напряжений под штампом 
ты а а г<1. 


Используя технику интегральных преобразований, поставленную задачу можно свести к решению 
следующего интегрального уравнения: 


1 со 
[т(р)р[Е (и).7 (игл) Л (ир^" ) диар = ^6 (0)гЕ, г <1, (1) 
0 0 
гдел =Н/а — безразмерный геометрический параметр задачи, характеризующий относительную 
толщину покрытия; .7, (х) — функция Бесселя 1-го рода; (и) — трансформанта ядра инте- 
грального уравнения, которая строится численно [6]. 
С использованием двухстороннего асимптотического метода решения парных интеграль- 
ных уравнений [7] можно построить приближённое решение задачи в аналитическом виде [6]: 


д 5(А^”) ) А (АХ), 


(г) =—=6 (0). 1, (0 нс я и ее т (2) 
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где постоянные С; определяются из системы линейных алгебраических уравнений: 


А В ыыы а т =0, К=1,2,...,М№; 


зай В? — А? [, (0) В? А? 
[м(и) — аппроксимация трансформанты ядра интегрального уравнения (1) выражением: 
м +” 
(и) = уве А,В, ЕС. (3) 


Решение (2) является асимптотически точным при Л -> 0 или Л -> 0 [7]. Точность решения для 
произвольного значения Л зависит от точности аппроксимации трансформанты ядра функцией 
(3). Алгоритм построения аппроксимаций высокой точности и связь между погрешностями реше- 
ния и аппроксимации описаны в работе [8]. 

Далее будем полагать, что функцию #(2), определяющую изменение модуля сдвига в по- 
крытии, описывает один из следующих законов: 











1) (2) = в не со$ (2п1^2), (= ю- А оз (2пк2), 
0 0 
_ |2 Е, 2 
35 (2) = 11 е[-н,0]/2.' 9* (2) 11 хе[-Н,0]/2,' 
г 47 +1, 4] +3 
тЫ 


В покрытиях 1—4 мягкие «слои» чередуются с жёсткими (рис. 1). Упругие свойства в покрытиях 1 
и 2 меняются непрерывно, в покрытиях 3 и 4 слои имеют чёткую границу. Параметр К соответ- 
ствует количеству периодов отрезка синусоиды, описывающей покрытия 1 и 2. Покрытия 3 и 4 
представляют собой пакет из (2А+ 1) слоя. В покрытиях 1 и 3 (2 и 4) модуль сдвига больше 
(меньше) либо равен модулю сдвига подложки на всём интервале 2 е [= ‚0 | 


О |1 и 1 1 





покрытия 1,3 ((5) покрытия 2, 4 
Рис. 1. Графики изменения модуля сдвига по глубине. Сплошные линии соответствуют покрытиям 1 и 2, 
пунктирные — покрытиям 3 и 4 


Показатель неоднородности К, = соп5{ >1 характеризует отношение модуля сдвига на 


поверхности Г к модулю сдвига подложки. Ограничимся рассмотрением случая Ё, = 7/2 , что соот- 


ветствует, например, сочетанию мягкого (алюминий, серебро, медь, свинец и т. д.) и жёсткого 
(железо, сталь, палладий, молибден и т. д.) металлов. 
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Анализ трансформант ядер интегрального уравнения. На рис. 2 изображены трансформан- 

ты ядер для законов 1—4 при К = 1, 10, 50. Из графика видна динамика изменения трансформан- 
ты ядра при увеличении числа слоёв. 

Пи) 

1.8 







о 


10° 0” 0” 10° 10. 10° 10° 10° 


Рис. 2. Трансформанты ядер для законов 1—4 при К = 1, 10, 50. Сплошные линии соответствуют законам 1 и 2, 
пунктирные — законам Зи 4 


Отметим свойства, общие для слоистых и непрерывно-неоднородных (функционально- 
градиентных, ФГ) покрытии: 
1) трансформанты ядер при всех значениях К имеют одну точку экстремума, обозначим её . Для 
законов 2 и 4 — это точка максимума, для законов 1 и 3 — минимума; 
2) при увеличении числа слоёв (параметр А) точка экстремума ш сдвигается вправо по оси и! 
3) значение экстремума трансформанты ядра [(ш) удаляется от 1 (для законов 2 и 4 увеличива- 
ется, для законов 1 и 3 — уменьшается) при увеличении К. При этом значения экстремума при 
К = 10и К= 50 близки, т. е. можно предположить, что они сходятся к некоторому пределу; 
4) при увеличении К расширяется диапазон значений и, в которых значение трансформанты су- 
щественно отличается от 1. 

Можно заметить ряд качественных отличий трансформант, соответствующих слоистым и 
ФГ покрытиям: 
1) при и < 0,5 трансформанты для законов 1 и 2 практически не изменяются с увеличением па- 
раметра КА, в то время как графики трансформант ядер для законов 3 и 4 сдвигаются вправо вдоль 
оси и, приближаясь к некоторому пределу; 
2) экстремальные значения трансформант (/(\%)) для слоистых покрытий, больше, чем для ФГ; 
3) для слоистых покрытий [(и) стремится к 1 при и -> со быстрее, чем для ФГ. 
Анализ контактных напряжений. Таблица 1 содержит величины погрешности аппроксимации 
трансформант ядер выражениями (3) для законов 1—4, где погрешность определяется формулой: 


А, = тах|[, (и)/Е(и)-1|.100 %. 


Таблица 1 
Погрешность аппроксимации трансформант ды для законов 1—4 
ВИННИ КВА: ВЯ ЗЫ СЕ ЗИ НЕ АЕ ИНН 98 


0, 16 % 0, 20 % 0, 34 % 


8% [03% | 
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Из значений таблицы видно, что построены аппроксимации высокой точности для всех 
рассматриваемых законов и значений параметра К, следовательно, погрешность построенных 
решений мала. 


Введём величину: Т,„ (Л,Г) =т(Л,Г)/Тьм (Г), где т(Л,г) — контактные напряжения под 
штампом для покрытий 1—4; Т,„(Г) — контактные напряжения, возникающие в однородном 


полупространстве с модулем сдвига, равным значению модуля сдвига на поверхности покрытий 
1—4, то есть С = (0). 

Величина т. (Л, г) наглядно демонстрирует отличие контактных напряжений, возникаю- 
щих на поверхности неоднородного и однородного полупространств, и, в отличие от т(Л, г), не 


имеет особенности при г = 1. Из физического смысла поставленной задачи следует, что для по- 
крытий 1—4 выполнено т,. (Л, г) > 1, при Л — 0 или Л -> оо. 


На рис. 3 изображена зависимость величины т, (Л, г) для покрытий 1, 2 при К = 1, 2, 10 


и покрытий 3, 4 при К = 10 от параметра Л в точке г = 0,5. Численный эксперимент показал, что 
графики для любого фиксированного Ге | 0..0, 99 | качественно не отличаются от приведённых 


на рис. 3. 
Те] ("| _ 


2.25 





законы 1, 3 
2 


Вэ" 


0.05 0.1 | 10 100 — 400 


Рис. 3. Зависимость относительных контактных напряжений от параметра Л при г = 0,5 


Отметим свойства контактных напряжений для покрытий 1—4, вытекающие из рис. 3: 


1) величины т,„ (Л,г) и т(Л,г) при фиксированном г имеет одну точку экстремума (обозначим 


ге! 
её Ло); 

2) при увеличении К точка экстремума Л сдвигается вправо по оси Л. Это означает, что макси- 
мальное отличие контактных напряжений для неоднородных покрытий 1—4 и контактных напря- 
жений, возникающих в однородном полупространстве, при увеличении К достигается для всё бо- 
лее толстых покрытий; 

3) при увеличении К значение экстремума функции т, (Л,„г) удаляется от 1 (для закона 2 уве- 


личивается, для закона 1 — уменьшается); 
4) при увеличении параметра К расширяется диапазон Л, в котором значение т, (Л, г) суще- 


ственно отличается от 1 (диапазон Л, где распределение контактных напряжений существенно 
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отличается от случая однородной среды). Другими словами, неоднородность оказывает суще- 
ственное влияние на всё более толстые покрытия. 

Свойства 1—4 являются общими для слоистых и ФГ покрытий, ниже приведены их основ- 
ные отличия: 
5) экстремальные значения контактных напряжений для слоистых покрытий меньше, чем для ФГ; 
6) кривые зависимости контактных напряжений от параметра Л для слоистых и ФГ покрытий по- 
хожи, однако, численные значения при фиксированных Л могут существенно отличаться. Напри- 
мер, при К = 10 иЛ = 18 отличие контактных напряжений для покрытий 1 и 3 достигает 25 %. 


На рисунке 4 изображены графики зависимости величины т, (Л,г)от координаты г для 


покрытия 1, А = 1 при «малых» значениях параметра Л (тонкие покрытия), «средних» Л (покрытия 
сравнимые с размерами штампа) и «больших» Л (толстые покрытия). 





ге] 
й 
А, = 0.631 
1.75 
— =2.24 
1.25 
А,= 5.25 
^,= 18.2 
0.5 0.75 0.99 0.5 0.75 0.99 0.75 0.99 
а) «большие» ^ 6) «средние» ^ в) «малые» А, 


Рис. 4. Относительные контактные напряжения для покрытия 1 при К = 1 


Контактные напряжения, возникающие на поверхности покрытия 1 при К = 1, обладают 
свойствами: 
7) контактные напряжения для «больших» значений параметра Л (Л > 1/2) монотонно убывают 
при г-› 1 по сравнению со случаем однородного полупространства. Для Л > 8 разница между 
т(Л,г) ит,» (Л,Г) менее 8 %, для Л > 25 — менее 1 %; 


8) контактные напряжения для «малых» и «средних» значений параметра Л (Л < 1/2) сначала 
возрастают, а затем убывают при г-+ 1 по сравнению со случаем однородного полупространства; 
9) из рисунка 4 видно, что при уменьшении Л сужаются области возрастания и убывания относи- 
тельных контактных напряжений по г (см. предыдущее свойство). Контактные напряжения для 
Л < 0,01 практически не отличаются от случая однородного полупространства везде, кроме 
окрестности точки г = 1, где наблюдается резкое возрастание, а затем убывание величины отно- 
сительных контактных напряжении. 

Следует отметить, что понятия «больших», «средних» и «малых» значений параметра Л 
условны, при изменении значения параметра К или рассмотрении покрытий отличных от 1—4 их 
числовые диапазоны могут существенно измениться. 

Свойства, аналогичные 7—9, имеют место и для покрытий 2—4 и всех значений параметра 
К (изменятся только количественные значения). 

В работе произведён анализ контактных напряжений, возникающих на поверхности покры- 
тий. Эта величина представляет наибольший интерес для теоретического исследования, т. к. зная 
её, легко можно рассчитать всё напряжённо-деформированное состояние полупространства [9]. 
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ТОВ$ТОМ ОЕ ЕГАЗТТС НАЕЕ-5РАСЕ МЛТН МОЕТТЕАУЕВЕО СОАТТМС ОЕ 
РЕВТООТС $ТВОСТОВЕ” 


А. 5. МУазПуем, Е. \У. Задаутп 
(Роп Чае Тесптса! Упмег$Ку), 
М. Е. УазПуема 

(Вооу Зае Тгапзрой Упмег$Ку) 


Ап ахсуттетс ргоМет о! {те ев5йсйу {Пеогу оп {огзоп оЁ ап е/азНс ва!-зрасе ит твотодепеоиз т аерт 
соайта о! репоас ягисиге Бу а поп-аеюгта Ме сисшаг рипсп 1/5 сопэ!аегед. Тре соайпд 1/5 гергезетеа Ву а $5ЕЁ 
ог аЁегпаНпа 50 апа Вага (ауегз. Тре Боипаапех о! {Пе Гауег$ сап Бе сеаг (р/есеил5е-сопапЕ ,атайоп о! пе 
5пеаг тоди/из) ог зто (сопйпиоизу твотодепеои$ ог Гипсйопа/у-дгааеа соайпсбс). ТВе арргохитаЕ апа/уй- 
са/ оивоп ю пе рго Мет оЁ ШАВ ассигасу 15 сопзёгисед и/Р/СВ [5 ейесйуе Гог апу соайпд и/сКпесх. Тре еТесЕ о! 
те питрег ог Гауег, апа о! {Пе соаНпа твотодепейу Ёуре оп Те Кете! гап5огт о! {пе тедга! едиайоп апа 
оп Пе сотасЕ $тез5ез @триноп ипаег те рипсй 1/5 пуезйдаеа. 1 15 5поип [ВаЁ иЁЛ себат иаиез оЁ те ргоВ- 
[ет рагатеегс, Пе зибзапйа!/ ОЙТегепсе Бейиееп те гези!Ёз Гог {пе сопбпиоиз апа р/есеилсе соп$апЕ Ёуре ог! 
пе уапабоп о! {те зВеаг тоди!из т аерт 15 оМатеад. 

Кеуигогаб: (огх/оп, ти[Пауегеа соайпдб, Гауегей сотро5[ес. 


" ТНе гезеагсн 15 допе мИЛ {Не зиррок гот {Не Миху оЁ ЕдисаНЧоп апа 5аепсе оЁ {Пе Виз$ап Еедега{оп (Зка{е Сопгасе 
Мо. 11.519.11.3015, Адгеетеп* №. 14.В37.21.1131) апа ВЕЕТ (дгагм{5 №05. 12-07-00639-а, 13-07-00954-а). 
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